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1. Introduction 

General relativity and quantum mechanics together imply that space-time structure at the 
Planck scale is not described by conventional notions of geometry. This was pointed out by many. 
In particular see Doplicher etal., [P as well as Werner Nahm (see [§]). These come about due to 
appearance of null horizons. Interestingly extra dimensional spacetimes with higher dimensional 
'Planck scale' 0(Tev) will require the extra dimensional space description to be fuzzy. Surpris- 
ingly the difficulties in defining geometry at infinitesimal distances were anticipated much earliar: 
'it seems that empirical notions on which the metrical determinations of space are founded, the 
notion of a solid body and a ray of light cease to be valid for the infinitely small. We are therefore 
quite at liberty to suppose that the metric relations of space in the infinitely small do not conform 
to hypotheses of geometry; and we ought in fact to suppose it, if we can thereby obtain a simpler 
explanation of phenomena ': Riemann, ^ 

Field theories on non-commutative geometries are inherently non-local leading to mixing of 
the infrared and ultraviolet scales. This, in turn, is responsible for new ground states with spa- 
tially varying condensates. Many non-perturbative studies have established that non-commutative 
spaces, such as the Groenewold-Moyal plane and fuzzy spheres, allow for the formation of stable 
non-uniform condensates as ground states. Exploring such implications of the non-local nature of 
field theories is very important in many areas of quantum physics. [Q, ||] 

Since different phases are intimately connected with spontaneous symmetry breaking (SSB), 
the role of symmetries in noncommutative geometries themselves is subtle. This issue is important 
in 2D because the Coleman-Mermin-Wagner (CMW) theorem states that there can be no SSB of 
continuous symmetry on 2-dimensional commutative spaces. There is no obvious generalisation of 
the CMW theorem for non-commutative spaces, since the theorem relies strongly on the locality of 
interactions. Non-commutative spaces admit non-uniform solutions (in the mean field) and one can 
ask the question what happens to the stability of these configurations. Non-uniform condensates 
naturally have an infra-red cut-off for the fluctuations. This cut-off softens the otherwise divergent 
contributions of the Goldstone modes[|, 0, |l0[ [TT], [| |, 0]. 



There have been various attempts to study gravity theories within the noncommutative frame- 
work [ f[5[ |l6| ]. This has led to a Hopf algebraic description of noncommutative black holes [17, 18] 



and FRW cosmologies [[T9|]. A large class of such black hole solutions, including the noncommuta- 
tive BTZ [^0, 21] and Kerr black holes, exhibits an universal feature where the Hopf algebra is de- 



scribed by a noncommutative cylinder [E2J, which belongs to the general class of the fC-Minkowski 



algebras [23, 24, 25, 26]. we shall take the noncommutative cylinder and the associated algebra as 
a model for noncommutative black holes. 

The study of quantum field theories in the background of black holes has led to the discovery 
of interesting features associated with the underlying geometry, such as the Hawking radiation and 
black hole entropy. In the noncommutative case, the black hole geometry is replaced with the al- 
gebra defined by the noncommutative cylinder. In order to probe the features of a noncommutative 
black hole, it is useful to analyze the behaviour of a quantum fields coupled to the noncommuta- 
tive cylinder algebra. Scalar field theories have been extensively studied on fC-Minkowski spaces 



[27, Eq, E9L 30, B3], which has led to twisted statistics and deformed oscillator algebra for the 



quantum field [29, 30, 31]. Theories on the noncommutative cylinder lead to quantization of the 
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time operator [33, 34]. See Madore [35] or Balachandran et al [[?!]] for an introduction to fuzzy 
geometries. 

In this paper we discuss examples of fuzzy geometries in Sec 2. In Sec 3 we consider QFT's on 
such geometries. Following this in Sec 4 we take up aspects of numerical simulations of flutuations 
of fields on fuzzy geometries and present our results. Lastly we conclude with discussions on 
implications of our results in Sec 5. 



2. Examples of fuzzy geometry 

It has been well known that representations of SU (2) Lie-algebra provide a basis for the study 
of functions on a fuzzy spheres. This can be understood by the quantisation of coadjoint orbits 
SU (2)/t/(l). The generators of the Lie algebra X; satisfy: 

[Xi,Xj] = ie ijk X k , £x 2 = Rl - t 2 - 1 ) 

9 1 

The fact 5 Z = CP 1 is a coadjoint orbit is useful in quantising this space. This can be extended to 
CP 2 = SU(3)/SU(2) ® £7(1) and any CP". The fuzzy torus is defined by the two generators U, V 
satisfying U V = e' e V U. Finite dimensional representations can easily be constructed for this 
algebra for rational 6. We will explain certain non standard ones. 

2.1 Higgs algebra 

The Higgs algebra is defined by = aZ + j6Z 3 , [X±,Z] = ^fX + . It can be easily 

checked that the Casimir for this algebra is given by: 

V = ±[{X+,X-} + g{Z)+g(Z-l)], (2.2) 

where g(Z) is 

g(Z)=C () + |z(Z+l) + |z 2 (Z+l) 2 . (2.3) 

For Co = ji 2 , (X = —2{2pL + 1), and j8 = 4 The Casimir reduces to the expression: X 2 + Y 2 + 
(Z 2 — ji) 2 . Equating the Casimir to 1 and plotting the function for different values of \i we see 
interesting topology change [p7[|. Such changes in the topology were first considered by Arnlind 




Figure 1: Surface plots depicting the change in topology for jj. = 0, 1, 1/2. 



etal [38]. 
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2.2 Fuzzy cylinder 

The NC cylinder is defined by the relation 

[Z, J* ] = ae^ (2.4) 

where Z is hermitian and e'^ 1 is unitary. Since we are interested in simulations, we have to discretise 
the above NC cylinder. 

In the rest of this paper, we will work with a = 1 without loss of generality, since the simple 
scaling Z — > Z/a can scale a away in the commutation relation (2.4). 

For this purpose consider the spin 7 irreducible representation (IRR) of the SU (2) Lie algebra, 
given by 

= 2X 3 , [X ± ,X 3 ] = T X±. (2.5) 

Since the operator Z generates rotations around the axis of the cylinder, it can be identified with 
X 3 . But when we use the finite dimensional representations of SU (2) we cannot implement fl2.4| ) 
with unitarity for e"^ . For this purpose, we decompose X + as product of a unitary and a Hermitian 
operator as given by 

X+ = R (2.6) 

Here R is necessarily singular and can not be inverted. However a partial inverse R can be found 
such that RR = P, the projector such that l—P projects to the kernel of 7?. Thus we get 

[Z,e^]P = e 1 * P (2.7) 

To find a representation for R and e'® , we can look at 

R 2 =X_X + = L 2 -L 3 (L 3 + 1) (2.8) 

which commutes with Z = L3. Remember that in the usual representation of angular momentum, 
L 3 \l,m>= m\l,m > with \m\ < I. Shifting the indices from to 21 + 1 = 27, we have j = m+l + 1, 
leading to 

X_X + \l,m> = [l(l+l)-m(m+l)]\l,m>=[(l+l/2) 2 -(m+l/2) 2 ]\l,m> 
X_X + \J, j > = {J 2 -(J- j) 2 } \J, j >= j(2J - j) \J, j > 

There is only one hermitian positive solution to this equation which takes the form 



R tJ = y/i(2J-i)8ij (2.9) 

which is diagonal as expected, and whose null space is along the top state \J,2J >. As a result, 
P = 1 - \J,2J ><J,2J\, and 

27-1 

R= Y, [i(2J-i)]~ l/z \J,i><J,i\ 
1=1 

where the sum stops at i = 2J — 1 so that there is a zero in the last position on the diagonal. It is 
now possible to deduce the first 27—1 lines of the unitary matrix e"^ from (2.6): 

>j* \JJ>=\J,j + l>, j<2J. 



X + \JJ> = y/j(2J-j)\J,j+l> 
X + \JJ > = t» R\JJ >= y/j(2J-j)J* \JJ > 
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Eq. (2.6) yields no equation for the last column which is instead determined from its unitarity. The 
columns 1 , • • • , 2J — 1 of e"^ , given by | J, j + 1 > , form an orthonormal set, as expected for a unitary 
matrix. Then the last column will be a vector orthogonal to all these vectors and thus can only be 
proportional to \J, 1 >. After normalisation that still leaves a U(l) freedom so that 

(e^ij = Sij +1 + 8 a 8 jaJ (2.10) 

where j3 can be any real number. For j8 = 0, e"^ is just a circular permutation of length 2J. 

2.3 Noncommutative BTZ blackhole 

We briefly summarize the essential features of a noncommutative black hole which is useful 
for our analysis. In the commutative case, a non-extremal BTZ black hole is described in terms of 



the coordinates (r,(j>,t) and is given by the metric [EQ, 21] 



ds 2 



Af- 



r~ 

I 2 



4r 2 



dt 



M 



4r 2 



dr 2 + r 2 



dtp 



2r 2 



-dt 



(2.11) 



where < r < oo ; — oo <£<oo 5 O<0<27T,M and J are respectively the mass and spin of the 
black hole, and A = — 1 jt is the cosmological constant. In the non-extremal case, the two distinct 
horizons r± are given by 



Ml 2 



1± 



1- ^ 



J 

Ml 



(2.12) 



An alternative way to obtain the geometry of the BTZ black hole is to quotient the manifold 
AdSi or SL(2,R) by a discrete subgroup of its isometry. The noncommutative BTZ black hole is 



then obtained by a deformation of AdS^ or SL(2,R) which respects the quotienting [17]. In the 
noncommutative theory, the coordinates r, <p and t are replaced by the corresponding operators r, 
and Z respectively, that satisfy the algebra 



ae" 



f,t\ 

2 







(2.13) 



where the constant a is proportional to £ /(r + — r_). We shall henceforth refer to ( |2.13| ) as the 



noncommutative cylinder algebra. Furthermore, i denoting the operator corresponding the the axis 
of the cylinder, it will be therefore identified as the operator Z in the following sections. 



It may be noted that the operator r is in the center of the algebra ( [2.1 3\ >. In addition, it can 
be shown easily that e ~ 2m / a belongs to the center of (2.13) as well. Hence, in any irreducible 
representation of (2.13), the element e ^ 27tl '/ a j s proportional to the identity, 

e -2nii/a =e iy^ (214) 

where y € R mod (2n). Eqn. ( 2.14J ) implies that in any irreducible representation of (2.13), the 



spectrum of the time operator t, or Z, is quantized y22[ , p3| , |34| ] and is given by 

spec t 



na 



2k 



(2.15) 



In what follows we shall set y = without loss of generality. 



The noncommutative cylinder algebra (2.13) belongs to a special class of K - Minkowski 



algebra and it appears in the description of noncommutative Kerr black holes [18] and FRW cos- 
mologies JT^]. We shall henceforth consider (2.13) as a prototype of the noncommutative black 
hole. 
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3. Scalar fields on fuzzy spheres/cylinder 

We shall now present our analytic and numerical analysis of scalar fields on different fuzzy 
geometries. 

Let <I> be a scalar field on a fuzzy sphere defined by spin j = (N — l)/2 representation. It is 
given by a.N x N matrix. We consider the action given by: 

S = — Tr {^MLi,®]} + R 2 {r<S> 2 + A<J> 4 } (3.1) 

It is easy to see the ground states are characterised by 

<I> = 0, and <£/0, but Tr® = (3.2) 

which corresepond to the uniform and nonuniform or stripe phases. We can obtain the continuum 
limit by taking N — > °°. The planar limit is obtained by taking R — > oo. One gets commutative 
planar or noncommutative planar (Moyal) limit depending on R 2 /N — > oo or finite. If we have 
a complex scalar field <I> then global U(l) symmetry can be broken contrary to the expectation 
from Coleman-Mermin- Wagner theorem in the NC limit. This is due to the nonlocality of NC 
geometries. This has been shown through simulations in [|TT]] . 

3.1 Topological aspects on fuzzy spheres 

But our interest in this work is to consider the topological aspects of nonlinear fields on fuzzy 
geometries. For this, we consider three hermitian scalar fields <!>,• with global 0(3) symmetry. The 
most general action upto quartic interactions takes the form, 



4k 



£|[Lf,*] \ 2 + R 2 {r\®\ 2 + iPe ijk ^j® k + X(\®\ 2 ) 2 +^\ [<&;,<&;] | 2 ) 



(3.3) 



In the mean-field the above theory admits a uniform condensate for r < 0. However the fluctu- 
ations of the Goldstone mode render this solution unstable. Apart from the uniform condensate 
the above model admits many meta-stable solutions. To simplify our arguments we consider the 
case j8 = =0: We are interested on those solutions which are stable due to topological 

obstructions. For example, 

®i = a L h with a = (3.4) 

The analog of this configuration in continuum space is the hedgehog configuration where the 0(3) 
spin vector on the sphere is pointing radially outward. The spin is parallel to the position vector 
on the sphere. This configuration is topologically stable as it cannot be smoothly deformed to a 
uniform one. Similarly the above configuration cannot be smoothly deformed to <I> = / which is 
also a solution. The above configuration corresponds to a winding number one map from the phys- 
ical space S F to the vacuum manifold which is also S F . All topologically stable configurations in 
the continuum limit, can be characterised by the second homotopy group ). For a discussion 
on topological classification of the maps S F — > Sj? see [12]. To study the net effect of topologi- 



cal nature of the background configuration and non-locality on fluctuations, we consider only the 
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winding number one configuration which is given in Eq.(3.4). As mentioned our plan is to compute 



these fluctuations numerically. Even before computing the fluctuations one can make some general 



remarks about the behaviour of the fluctuations Ql 1|]. The effect of nonlocality basically provides 



a non-zero mass 0(a/N) to the Goldstone mode fluctuations. This puts an infrared cut-off for the 



fluctuations. From our previous study [ |1 lfl it seems that this mass/cutoff is mode dependent, as 
only higher modes of the condensate survived the fluctuations. As we will see from our results the 
combined effect of topology and nonlocality, the infrared cut-off drastically reduce the contribution 
of the fluctuations. Before presenting details of numerics let us also consider scalar field on fuzzy 
cylinder. 

3.2 The Action on the fuzzy cylinder 

Define TrO = Tr{POP) where P is the projection operator defined in [Eq: No] This trace 
TV is equivalent to integrating over the whole cylinder in the continuum limit. We also need the 
derivatives d§ and dz- They are: 

= [ZM (3.5) 
d z <D = e~* (3.6) 

Then, apart from /-dependent normalisation factors, the action can be chosen as: 

S = TV (| [Z,<D] | 2 + | e-'V*,*] I 2 + V(*)) (3.7) 

where V(<t>) is the potential which can be taken to be of the form, 

y(<I>) =n<t> 2 + c <t> 4 (3.8) 

for a hermitian field <J>. 

This action has a problem of instability coming from 7V(<!> 4 ) = Tr((P®)® 2 ($>P)) which 
cannot contain any quartic (nor cubic) term for the variable <t>2J2j- This makes the theory unstable 
with respect to this variable. The simplest cure is to insist that this term is not a degree of freedom 
of the theory and constrain it to be zero. To keep the set of fields an algebra, we set to zero the last 
row <i>2Ji and column <J>,2/ of the field. As a result the hermitian field <1> now only has (27 — l) 2 
degrees of freedom, and <I> = P<t>P. 

3.3 Dimensional reduction 

With this new choice of the field the action becomes 

S = TV( P | [Z,P<i>P] \ 2 P + P\ e-*[<P \P&P] | 2 P+ V(<D)) 
= TV (| [PZP,<t>] | 2 +| [Pe^P,®] | 2 + V(<t>)) (3.9) 

which can be rewritten simply as the action for a hermitian matrix in a (27 — 1) X (27 — 1 ) matrix 
algebra of reduced dimension: 

S = Trji (\ [Z,4>] | 2 +| | 2 + V(<D)) (3.10) 
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where f = J — 1/2 is the reduced angular momentum, while e l § and Z are the matrices obtained 
from e"^ and Z by removing the last line and column. For e'$ , it is equivalent to setting e'^ — > 
in its 2/ -dimensional expression ( 2.1C ). As for Z, it is therefore the 2J' x 2J' diagonal matrix 
obtained from Z by removing its top eigenvalue J — 1/2: 



Diag(-/+ 1/2, —J + 3/2, 3/2) = Diag(-/, -/+ 1, 1) 



■Zij = (—/ — 1 + i)8ij. 



(3.11) 



Note that Z and Z are defined by their commutation relation ( |2.4[ ) and only appear in the action 
through dip as a commutator. As a result, they are only defined up to a translation by a matrix 
proportionnal to the unit, and thus Z = Diag(l, • • • ,2J') is another possible choice. 
Although is not unitary, the equation 

< m\[Z,e"l']\m > = (m —m)8 m ^ m+ \ = 8 m i >m+ \ =< m'\e'^\m > if m < 2f 
= 0=<m'\[Z,7}\2J' > ifm = 2/ 



shows that Z and e"r do satisfy the commutation relation ( J2.4| ) for a = 1 

The cylinder is also parametrised by its radius r. According to ( 2.13 ), r commutes with both Z 
and e'^ . It can therefore be considered as a pure number in the non-commutative cylinder algebra. 

The radius will appear as a simple scaling in the action. The volume of the cylinder 7Y(1) 
depends linearly on r, so the action should have an overall scale of r. The derivative along the axis 
dz does not scale with r, whereas the angular derivative d$ scales like 1 jr. 

As a result, the action on a fuzzy cylinder of radius r is given by: 

1 



(3.12) 



rTr[-,\ [Z,<D] \ l + | ^[^,<I>] \ l + V{$>) 



r Try ( i| [Z,<J>] | 2 +| [<?'"*, <D] | 2 + V(4>) 



3.4 Spectrum of the Laplacian 

The spectrum of the Laplacian can be obtained from: 



1 



t 



if 2 0= [Z,[Z,<I>]] + -([<?'"* ,[«?'■*,<&]] + ,0]]) 
= ^ z 2 0+^(^_^ + + ^ + ^_<I>) 

= ^ z 2 <J> + ^_^ + + i[[^,^ t ],<I>] =^ z 2 <I> + ^_^;<I>+^[Diag(-l,0, 
The eigenmatrix equation for the Laplacian reads simply 

££ 2 Q m (d) = X& m (d) = ® m (Xd) <^M m d = Xd 

where 

2J'-m 

0> = ® m (d)+®l l (d), v/ith® m (d) = di\i><i + m\,0<m<2J'-l 



(3.13) 
,0,1), *] 



(3.14) 
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for the eigenmatrices, and the vector d = (d/)i<;<z/'-m is the unknown to be determined by the 
eigenmatrix equations. 

For the Laplacian on a cylinder of radius r, we get 



/ 1 -1 
-1 2 -1 



\ 



(0) 



V 



(0) 

2 -1 
-1 1 



/ 3/2 -1 
-1 2 -1 



m 



+ 



\ 



(0) 



(0) 



V 



2 -1 
-1 3/2/ 



, m^O (3.15) 



These matrices are similar to the ones obtained for the Laplacian on a one-dimensional lattice and 
can actually be diagonalised without much difficulty, taking good care to remember that M m is a 
matrix of dimension 2J' — m. 



Spectrum of Mo The eigenvalues of Mo are: 

Xq = 4sin 2 (*7r/4/ , ), < k < 2f. 
Spectrum of M m , m / Reparametrizing the eigenvalues as 

A =2-2cos(0) =4sin 2 (0/2), 



(3.16) 



(3.17) 



We can solve for eigenvalues for any N numerically. However for large matrices N ^> 1, it is 
possible to find approximate ones: 

• For 6 < %, or k < N, tan((7T - 0)/2) ~ 2/0 > 1. Therefore, Nd = k% + n/2 - p k with 
Pk <C 1. The equation then becomes: 

^" 3(^ + ,/2-p fc) ^ P ^2W^ + 1/2) - 

• For ~ ft, or £ ~ N, tan((ft - 0) /2) ~ (ft - 0)/2 < 1, and therefore, NO = kn + p k with 
Pi <C 1. The equation then becomes: 

A/ft -k%- p k N -k 

Pk - 777 ^ P-t - "^r^TT 



6iV 



which is a small number, as expected, since k~ N. 



6N 



4. Numerical Simulations and Results 

Effects of the fluctuations beyond mean field are computed from the partition function, which 
in the path integral approach is given by, 

2? oc J D<t>e- S ^l (4.1) 
The standard numerical methods adopted for this integration are Monte Carlo simulations. 
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4.1 The numerical scheme: pseudo heatbath 

In the Monte Carlo algorithms, one generates an "almost" random sequence of <J> matrices 
by successively updating elements of <J> taking into account the measure and the exponential in 
the integral above. This sequence of <J> is then used as an ensemble for calculating averages of 
various observables. For a good ensemble the auto correlation between the configurations in the 
sequence must be really small. Though this auto correlation can be reduced by using some over 
relaxation programme [g], it is greatly reduced, however, when "heatbath/pseudo-heatbath" type 
of algorithms are used. This method is very much common in the non-perturbative study of <J> 4 
theories in conventional lattice simulations. It gives better sampling and is efficient at least for 
smaller A values. This is why we make use of "pseudo-heatbath" technique [|l0|,[ll|]. 

4.2 Topological stability and 0(3) model 

In our simulations, for each choice of parameters, we choose an initial configuration given by 



Eq. ( |3-4| ). Fluctuations around this configuration are then generated by the above updating method. 
Since this configuration is a variational solution to minimising the classical action, it will thermalise 
as we update/include the thermal fluctuations. Once the initial configuration is thermalised we 
compute the observable M. We make measurements after every 10 updates of the entire matrix. 
We also use over-relaxation to reduce the auto correlation of the configurations generated in the 
Monte-Carlo history. 



In a numerical simulation, the condensate will not maintain its exact form as in Eq. (3.4) along 
the Monte Carlo history. The configuration can evolve into different random SU (2) rotated config- 
urations of Eq.([3~4|) as we keep updating it. To overcome this, one needs to rotate the configuration 



at each step of the Monte Carlo history so that the configuration takes the form of Eq. ( |3.4| ). But 
this is a difficult and time consuming task. On the other hand one can have an observable made of 
<I>-s which is invariant under the SU (2) rotations, e.g basis independent. For this purpose we define 
the following observable, 

Aij = ^Tr(L^j),M = VAU (4.2) 



M projects out the / = 1 angular momentum mode. Note that the initial configuration in Eq. ( p.4| ) 
projects out only the 1=1 mode. Analysing the statistical behavior of M will give us a definite 
conclusion about the stability of the initial configuration. We mention here that Tr(£&f) is also an 

i 

SU (2) invariant. But the information on the amplitudes of different / modes gets lost in this form. 
Also comparatively the observable M may serve as an order parameter in the case of any phase 
transition of the hedgehog configuration to <J>, = at high temperatures. We mention here that 1 = 1 
is the lowest possible stable mode, as / = mode will be unstable. One can consider configurations 



with higher winding, instead of Eq.( |3.4[ ), however we expect them to be more stable than the / = 1 
condensate. This is because the infrared cut off will rise with higher winding configurations. 

For practical reasons, the size of the matrix N, in other words size of the resolution scale is 
finite. So there are usually finite volume (R,N) effects. So a non-vanishing condensate <!>,• does not 
mean there is SSB. One needs to define suitable observable dependent on <!>, which should scale 
with (R,N) appropriately in the thermodynamic limit (N — > °°,R — > °°) to conclude anything. Now 
there are two possible thermodynamic limits. If in the thermodynamic limit the ratio R 2 /N does 
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not vanish then the space is described by a non-commutative algebra. This limit is of interest to us, 
as we expect that the CMW theorem will hold good in the commutative thermodynamic limit. 
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Figure 2: Monte Carlo history for N=48 
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Figure 3: Histogram H(M) for N=48 



4.2.1 Results and discussions 

In our calculations we fix R 2 /N = 10, r = —8. For simplicity we take X\ = 0.25. With 
this choice of parameters we do our simulations for five different sizes of the <I> matrices, N = 
48,56,64,78,96. Fig.g gives a typical Monte-Carlo history of our simulations for N = 48. 

In FigJ^ M fluctuates around a value close to the initial value. Then M suddenly jumps to a 
small value and settles down. A histogram H(M) of M clearly shows two peaks M as seen in Fig.||. 
The peak on the left has large / = and small 1 = 1 component. The peak at higher value of M has 
large / = 1 component and small / = 0. This peak is close to the value of the initial configuration. 
So in this state fluctuations modify the initial configuration slightly and retain its topological nature. 
In our Monte-Carlo history we observed the / = 1 state decaying to / = state but not vice-versa. 
This implies that due to finite volume effects, the uniform condensate is more stable than our initial 
hedgehog configuration for this case of N = 48. 

To study the stability of the 1 = 1 configuration we considered both the commutative and non- 
commutative limit. For the commutative limit we fixed R and considered higher values of N. We 
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did not observe any change in the distribution of M in the 1 = 1 state. The average value, and the 
fluctuations of M remain almost the same as we go from N = 48 — > 64, as can be seen in Fig.|]. 
As for N = 4$ the 1 = 1 configuration also decays for N = 64. This result suggests that the 1 = 1 
topological configuration is not stable in the commutative continuum limit as expected. 




4000 8000 12000 16000 20000 

M-C history 



Figure 4: M-C history for fixed R 2 



1.2 
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Figure 5: M-C history for fixed R 2 /N 



There is a complete change in the behavior as we consider the non-commutative limit, i.e 
fixed R 2 /N as we increase N. Except for the lowest N = 48 the 1 = 1 state did not decay during 
the entire run for higher N. In Figj5| we show the Monte-Carlo history of N = 48,64,96. Unlike 
the commutative limit, the fluctuations of M decrease with N. In Fig|6|, we give the average value 
of M as a function of N. The average value of M increases slightly with N, with the variation 
decreasing with N. This suggests M will reach a finite value in the continuum limit. We also 
compute the fluctuations of M to see any possible scaling with the cut-off N. In Fig.0, we show 
X = (M 2 ) — (M) 2 in the / = 1 state. The solid curve represents a fit, f(N) ~ A^ a with a ~ —4.. 
This clearly suggest that the 1 = 1 state is stable in the N — > °° leading to spontaneous breaking of 
the 0(3) symmetry. 

We also mention here that one can start with an initial uniform 1 = configuration and consider 
fluctuations. We expect that the results be similar to that in ref.[ll |. In ref.[ll ] it was found that 
only the highest mode / = (N — l)/2 condenses. The fact that we find the 1 = 1 mode stable clearly 
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Figure 6: (M) vs N 
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Figure 7: % vs N 



shows that the topological nature of the initial configuration complements the effect of non-locality. 
These two effects drastically reduce the fluctuations. 



4.3 Fuzzy blackholes: NC cylinder 



The model defined by the action Eqs.( |3.8y3.12| ) which we want to simulate has three parameters 
(jU,c,r) plus the matrix size J. The goal is to explore the parameter space for various phases of <J>. 



The simulations are carried out using the "pseudo-heat bath" Monte-Carlo (MC) algorithm [10, 11] 
to reduce the auto-correlation along the MC history. 

The field should also be allowed to explore the whole phase space and not remain trapped in 
local minima. To this end, an over-relaxation method, first suggested in is also used. Let us 
introduce S<j>(<I>; ; -) the dependence of the action on the field entry <J>,y when the field takes the value 
<!>. It is a fourth degree polynomial. Therefore the equation S<j> (<£>,-,•) = S®(<&jj = a), which has an 
obvious solution <J> iy = a, can be factorised into a degree three polynomial which always admits 
at least one real solution. The overrelaxation method consists in replacing the field entry <J>, 7 = a 
by one of these real solutions, thereby moving the field in a different region of the phase space. A 
crosscheck is also used to verify that the field probability distribution of our Monte-Carlo runs are 



13 



Phase structures in fuzzy geometries 



T R Govindarajan 



consistent. Let us split the terms in the action according to their scalings 

S(0) = 5 2 (0) + S 4 (4>) with Si{X$) = X%{<S>). 
Then one can define a modified partition function 

Z(A) = y [d0]e" 5(A0) = j [ d( j)] e -^S2W-^s4i) (4 3) 



X~ N j[&y]e- s ^\ yf = X(j>, (4.4) 



where is the number of degrees of freedom in the field <p which appear in the integration. Evalu- 
ating 

dln(Z) 



= -2 < S 2 > -4 < S 4 > from 
= -N from (0) 



yields the check originally due to Denjoe O'Connor [p9Q. 

<S 2 >+2<S 4 >=N/2 (4.5) 
In all simulations, this identity ( |Q1 ) is always satisfied to better than 1% relative error. 

4.3.1 The phase structure 

The temperature (J) is regulated by varying the parameter /j.. 

• /I <C 1 corresponds to low temperatures when the fluctuations are small. In this case, the 



minimum of S gives the most probable configuration of the phase. In Eq. (3. 12), it is possible 
to minimise the action by minimising separately the kinetic term, so that <& « 1, and the 
potential term so that <I> = \J —\ij1c\, and this phase is therefore known as the uniform 
phase. 

At high temperatures, \i » 1, the thermal fluctuations lead the system to the disorder phase 
<I>~0. 

At intermediate temperatures, the competition between the action and the fluctuations give 
rise to new phases called the non-uniform or stripe phases. These new phases are specific to 
non-commutative spaces. Various numerical studies have confirmed the existence of these 
phases [^, [7|, ^, K], 11, |, 14, 13] on the fuzzy sphere. A non-commutative cylinder will also 



exhibit the non-uniform phases. However, due to the non-trivial topology of the cylinder 
(the first homotopy group being non-trivial), one can have a more complex phase structure 
described below. 

For example there can be stripes going around the cylinder, or parallel to its axis. These two 
phases can be distinguished by their overlap with the operators Z, e 1 ^ , and e 1 ^ respectively. Stripes 
going around the cylinder will have non-zero overlap with the operator Z. While a configuration of 
stripes along the axis will have overlap with e 1 ^ and e 1 ^ . We present our results in the following 
subsection. 
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4.3.2 Example numerical runs 

For a given choice of N = 1, c = 0.36, and r = 1 the simulations are done for various values 
of ju. The various phases discussed above can be characterised by the observables m u = Tr(<&), 
m z = 7>(<1>Z), m x = Tr(<&e"^). A finite m u with (m z ,m x ) ~ characterizes the uniform phase. 
On the other hand, (m u ,m x ) ~ with non-zero m z characterizes stripes going around the cylinder. 
Stripes along the cylinder charactersied by (m u ,m z ) ~ with non-zero m x . 

For jX = —35.1, the data of a run are shown on Fig.||, and, as expected, we observe the uniform 
phase. 

For n = —20.0, we observed the phase with stripes going around the cylinder. This is verified 
on the histogram of the observed values of m u ,m z plotted in Fig.|[ It is clear from the figure that 
the average value of m z is finite while the average value of m u is vanishingly small. 

Fig.[l0| shows the system in the disorder phase where m u ,m z ,m x all fluctuate around zero. We 
did not observe the phase with stripes going along the cylinder as a ground state for any choice 
of /I for r ~ 1 . One can expect to observe this state for very small r when the second term which 
suppresses this state is made subdominant. For a very small radius r = 0.01, c = 36. and /i = 
—3.6 x 10 3 , this phase appears as meta-stable in Fig.|IT[ This phase is stable w.r.t small fluctuations. 
Only large fluctuations, which occur less frequently, can destroy such a state. 
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Figure 8: m UlZ , x vs MC history N = 7,n 
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Figure 9: Histogram of H(m u ),H(m z ) :N = 7,11 = -20 
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Figure 10: Histogram of H(m u ),H(m z ) : N = 7 = 0.36 
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Figure 11: m u , z , x vs MC history jV = 7,£i = -3.6xlOV = 36,r=.01 



5. Conclusions 

In this analysis we have shown that topologically non-trivial configurations on the fuzzy sphere 
avoid the CMW theorem much more dramtically than the non-topological symmetry breaking [ftOj]. 
The mass gap or the infrared cut-off in this case is large enough to render the fluctuations of the 
Goldstone modes finite. On the other hand for non-topological condensates the Goldstone modes 
are large enough to destroy almost all the modes except the few highest modes. We have presented 



the simulations wherein the cubic Chern-Simons (CS) term is absent in the action Eq (p.3|). The 
Chern-Simons term allows topological solitons even when the quadratic mass term is positive upto 
some value. Interestingly with CS term the configuration 0, = a Li is preferred over symmetric 
solution. On the otherhand, it is not expected to alter the picture of topological stability of the 
solutions. This term plays an important role in the emergent geometry in NC fuzzy spaces [ p| , |42| ] . 
What we find here in the simulations is that even in the absence of CS term, emergent fuzzy spaces 
can be stable. The stabilty of higher dimensional fuzzy spaces like CPp are of significance in 



this context Q43[] - The implications of this stability for extra-dimensional fuzzy spaces will be 
considered later. 

We have also considered a finite dimensional representations of the noncommutative cylinder 
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algebra, which make it fuzzy. We study scalar field theory in the background of this algebra both 
analytically and using numerical simulations. 

In the numerical simulations of scalar field with a generic potential we find, as expected in 
noncommutative cylinder, novel stripe phases breaking translational symmetry. But they have some 
differences with the usual stripes on Moyal spacetimes. These are also stable due to topological 
features arising in this fuzzy geometry. 

It is well known that a large class of noncommutative black holes are described by a non- 
commutative cylinder algebra. The fuzzy cylinder algebra derived from it can therefore be used 
to define a fuzzy black hole. From general considerations [|l]] we know that such black holes can 
arise at the Planck scale. Our results provide a first glimpse about the phase structure of a quantum 
scalar field theory in the background of a fuzzy black hole at the Planck scale [44]. 

Acknowledgment This work was done as a part of the CEFIPRA/IFCPAR project 4004-1 en- 
titled Fuzzy Approach to Quantum Field Theory and Gravity. The authors gratefully acknowledge 
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